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Abstract 

Recent numerical work on the fate of plasma instabilities in weakly-coupled non-Abelian gauge 
theory has shown the development of a cascade of energy from long to short wavelengths. This 
cascade has a steady-state spectrum, analogous to the Kolmogorov spectrum for turbulence in 
hydrodynamics or for energy cascades in other systems. In this paper, we theoretically analyze 
processes responsible for this cascade and find a steady-state spectrum ~ k~ 2 , where is the 
phase-space density of particles with momentum k. The exponent —2 is consistent with results 
from numerical simulations. We also discuss implications of the emerging picture of instability 
development on the "bottom-up" thermalization scenario for (extremely high energy) heavy ion 
collisions, emphasizing fundamental questions that remain to be answered. 



I. INTRODUCTION 



It is important to understand theoretically the mechanisms by which quark-gluon plasmas 
can locally equilibrate in heavy ion collisions such as those at RHIC. Since the path-breaking 
work of Baier, Mueller, Schiff and Son p| on "bottom-up" thermalization, one modest the- 
oretical goal has been to understand the process of equilibration in the simplifying theoret- 
ical limit of arbitrarily high-energy collisions, where the running strong coupling a s can be 
treated as small. Typically, theorists expect to be able to solve weakly-coupled problems, 
but equilibration of weakly-coupled non-Abelian plasmas has proven to be a very rich and 
challenging problem. In particular, the original attempt by Baier et al. did not account for 
the physics of plasma instabilities, which are now believed to play a crucial role during some 
of the early stages of thermalization 0, 0| • 

Early in a heavy-ion collision, particles have an anisotropic distribution of momenta, as 
measured in local frames moving with the expanding plasma 0. Later, they scatter and 
equilibrate to locally isotropic distributions, giving rise to hydrodynamic behavior. Generi- 
cally, anisotropic distributions of particles produce collective plasma instabilities known as 
Weibel or filamentary instabilities 0, 0, 0,11, EiJ U|- These instabilities are associated 



with filamentation of particle currents and initially exponential growth of long- wavelength 
magnetic fields. Since these magnetic fields in turn deflect the particles in the plasma, and so 
perhaps drive the particles towards isotropization and equilibration, it has been imp ortant 
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pernaps drive tne particles towards isotropization and equilibration, it nas been impc 
to understand just how large these unstable magnetic fields grow 12|, H, HI IE la J 

In this paper, we will follow the nomenclature of Baier et al. and refer to the initial 
particles of the plasma {e.g. small x gluons in the saturation picture) as "hard particles," 
and we will refer to their momentum scale {e.g. the saturation scale) as Phard- We will 
refer to the momentum scale of magnetic plasma instabilities as m so f t . Once the system has 
expanded enough that the hard partons have perturbative density (nhard <C Phard/^ 2 ); one 
generically finds that m so f t <C Phard jH 03 • We will assume and exploit this hierarchy of 
scales in what follows. We will refer to any scale that is parametrically small compared to 

Phard as "Soft." 

Fig. shows a cartoon of the results of numerical studies in 3+1 dimensions of non- 
Abelian plasma instability growth 0, 0] in a non-expanding system. The energy in soft 
magnetic fields is plotted versus time for simulations initiated with a very tiny magnetic 
seed field. The soft magnetic field energy initially grows exponentially, but the behavior 
changes once the fields become large enough that their non-Abelian self-interactions become 
non-perturbative {B 2 ~ m^ oit /g 2 ). The energy growth then becomes linear with time. For 
a non-expanding system, this linear growth would continue until some parametrically long 
time later (not addressed by the simulations) when the soft fields eventually have a non- 
perturbatively large effect on the original, anisotropic hard particles that created them. 

Recent numerical work 0] indicates that the linear growth in energy is associated with 
a cascade of energy from the soft, perturbatively unstable modes toward the ultraviolet. 
A cartoon of the development of the spectrum is shown in Fig. |21 (a) The system starts 
with the initial, anisotropic hard particles of momentum ~ Phard plus some small initial 
fluctuations in softer modes, (b) Plasma instabilities initially lead to rapid, exponential 
growth of unstable soft modes of momenta ~ m so { t . (c) The unstable modes stop growing 
in amplitude once they become non-perturbatively large, but total soft energy continues to 
grow (linearly) as interactions feed power from m so ft into increasingly higher momentum 
modes. The spectrum of this cascade falls with some characteristic power k~ u , and the 
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FIG. 1: A qualitative depiction of soft magnetic field energy vs. time for the growth of non-Abelian 
magnetic instabilities. 
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FIG. 2: A cartoon description of the energy spectrum of excitations at various times during the 
development of non-Abelian plasma instabilities. 



upper momentum limit A(£) of the cascade increases with time. The goal of the present 
paper is to theoretically determine the exponent v by considering the processes responsible 
for this cascade. 

The basic plan of our analysis will be somewhat similar to investigations of thermalization 
during reheating after inflation in the early Universe ^1 20 , 21 1 . The details of our problem, 



and the result for the spectral index u, are different, however. 

Let fk be the average occupancy of each mode k. In analyzing the cascade, we shall take 

1 « /* « ~ 2 (1.1) 

9 

for m soft <C k <C A(i). In this range of /, the physics is perturbative (fk <C V^ 2 ) an -d can 
be equally well described by a kinetic theory of particles or by a classical theory of fields 
(fk 3> 1) [13 - 1 The perturbative nature of the cascade for k 3> m so{t follows from (i) the 
non-perturbative occupation number / ~ 1/g 2 at the instability scale m so f t and (ii) the 



This is a slight oversimplification. Perturbative occupation numbers fk <C 1/g 2 are necessary but not 
sufficient for a kinetic theory description of relativistic systems. For quantitatively accurate (or qualitative 
order-of-magnitude) applications of kinetic theory, one must check that deBroglie wavelengths and the 
duration of scattering events are <C (or <) mean free paths, for whatever processes dominate the kinetic 
description. 
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decrease of / as some power k u to be determined. This gives 

, ss 1 fm^y , . 

where the superscript "ss" stands for "steady-state." Perturbative cascades are known in 
the plasma physics literature as "weak plasma turbulence" We will restrict attention 
to the case where A <C g~ 2 ' v m B0 ft in order that (jl.2|) satisfy the other condition ^> 1. 
This classical field theory limit is implicit to the numerical simulations which discovered the 
cascade [3, 17, 3], since these simulations treated all soft fields as classical. 



In the next section, we warm up by discussing 2^2 scattering of particles at the UV 
end of the cascade. We show that this process alone would reproduce a result for the 
spectral index v in (jl.2j) known from analogous scalar theory applications (v = 5/3). In 
Sec. IIIII we show that there is a more important process: multiple scattering from the 
non-perturbative infrared background produces a steeper distribution, v = 2. This result is 
consistent with numerical simulations [l8[. Sec. II VI discusses the back reaction the cascade 
has on the infrared physics of the instability and poses some open questions concerning 
logarithmic effects in our analysis. In Sec. 13 we return to the question of how instabilities 
affect the bottom-up scenario and discuss what remains to be understood about instability 
development before this question can be answered. Finally, we offer our conclusions. 

Non-perturbative numerical studies of instability development in 3+1 dimensions have 
only studied hard particle distributions which are moderately, but not extremely, anisotropic 
TlHEi. In the case of moderate, 0(1) anisotropy, there is a single scale which charac- 
terizes m so f t , given parametrically by 0,11(1 11 

p xhard 

<*~<? r-*—, (1.3) 

where the integral is over the distribution of hard particle momenta p. For extremely 
anisotropic distributions, the situation is more complicated. In most of this paper, we 
will implicitly assume that there is only a single relevant soft scale. We will only discuss 
the unresolved issues surrounding extremely anisotropic distributions in our discussion of 
bottom-up thermalization in Sec. |3 



II. WARM-UP: HARD 2^2 SCATTERING 

The standard picture of driven turbulence is that something pumps energy into the 
cascade from the infrared (IR) end at a constant rate. In the absence of dissipation, this 
leads to linear growth in the total energy. At the same time, energy is flowing out the ultra- 
violet (UV) end by increasing A. The spectral index v of the cascade can be determined by 
balancing (i) the rate at which energy flows out the UV end of the cascade by scattering 
processes that extend the UV cascade by increasing A (e.g. scattering momentum A particles 
to produce momentum 2A particles) and (ii) the constant rate energy enters the IR end of 
the cascade. 

As a review and warm-up example, let us for the moment (incorrectly) assume that the 
dominant process for increasing the cut-off was large-angle 2— >2 scattering with all momenta 
roughly of order A. Fig. El shows an example of how two particles with momenta p\ and P2 
less than A can produce a particle with momentum greater than A in 2^2 scattering. Now 
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consider what the time scale would be for bringing f 2 \ up to its steady-state value (|1.2jl . 
For a statistically homogeneous system, the Boltzmann equation for 2— >2 scattering has the 
form 

4/(P4)= / l-Mp lP ^ P 3P4| 2 [/l/2(l + /3)(l + /4)-/ 3 /4(l + /l)(l + / 2 )] (2.1) 



for bosons such as the soft gluons that make up our cascade. In the / 3> 1 classical field 
limit of (jl.ip . this becomes 



d 
dt 



/(P4) = / |-M PlP2 ^ P3P4 r[/i/ 2 (/3 + A) - / 3 / 4 (/i + / 2 )]- (2.2) 



Even if f{p±) were initially small for p^ > A as in Fig. there is a term in this equation with 
/1/2/3 ~ /a- Take the square-amplitude |.M| 2 for large-angle scattering as order g 4 times 
normalizations and the energy-momentum conserving delta function. The characteristic rate 
T implied by (|2.2|) for /(2A) to grow to its steady state value /| A ~ / A s is then given by 

r/r ~ gVTfA, (2.3) 

where the factor of A on the right-hand side follows from dimensional analysis. So 

T ~ <? 4 (/ A s ) 2 A. (2.4) 

From our assumed power-law form (jl.2|) for / ss , this gives 

T oc A 1 ' 2 ". (2.5) 

Since the cutoff A has to change on the time scale of the age t of the system, we must have 
r ~ l/t. 

On the other hand, let us estimate the rate from the linear growth of the total energy 
density in the cascade. Using (jl.2|) . the total energy density is 

fty fc ^(^) V ocA 4 - v (2.6) 

(provided v < 4, which we will find). Taking e oc t, the corresponding rate is 

j oc A u ~\ (2.7) 

Comparing the exponent here to (|2.5|) . we would obtain v — 5/3, in agreement with previous 
analysis in the literature of certain cascades in field theory where large-angle 2— >2 scattering 
is the dominant process 
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2 In particular, this v agrees with Eq. (66) of Ref. [ill ]. 
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FIG. 3: Momenta in a P1P2 — > P3P4 process that creates a particle with momentum larger than 
A. 

III. THE DOMINANT PROCESS 

The more efficiently the system can scatter energy into the UV, the larger will be the 
exponent v. The goal is to identify the most efficient process. We propose that the dom- 
inant process is multiple scattering of quanta of momentum A off of the non-perturbative 
background field of wave-number m so f t . In this section, we analyze this process and show 
that it is more important than the large-angle 2— >2 scattering considered in the last section. 

Consider a relatively high-energy particle {k ~ A m so f t ) propagating through a soft, 
random, non-perturbatively large background field Abkgd having wave- numbers of order m so ft. 
The particle can scatter from the background as in Fig. 0] Non-perturbatively large means 
(Mbkgd ~ "^soft, and so the only parameters in the problem are k and m so ft. The probability 
of interaction is independent of k in the k 3> m so ft (Eikonal) limit. By dimensional analysis, 
the rate is therefore O(m so f t ) to change momentum by O(m so f t ). Each such scattering will 
have a perturbative effect on the particle, randomly changing its momentum k by m so f t <C k. 
N successive collisions will change k by order N l / 2 m so { t in a time of order N/m so f t . Taking 
N ~ (k/m so{t ) 2 ~ (A/m so f t ) 2 , we see that the time taken to change particle momenta from 
A to 2A at the UV end of the cascade would be 

*~(^V^, (3.1) 
\m soit J m soft 

corresponding to a net rate that depends on A as 

r A ^2A ~ - oc A" 2 . (3.2) 

If we equate this to the rate ()2.7|) determined by linear energy growth, we obtain the spectral 
index 

v = 2. (3.3) 

The fact that this v is larger than the v = 5/3 of Sec. |H] indicates that the large- angle 2^2 
processes of Fig. El are unimportant compared to multiple scattering by Fig. HJ 3 

We note in passing that interactions with a non-perturbative background field have also 
been important in the study of cascades during reheating of scalar field theories in the 



3 One can check this explicitly by comparing the rates for both v = 5/3 and v — 2: Eq. 1)2.4(1 gives 
"V > ft(?Ti S oft/A) 2ly ~ 1 for large-angle 2— >2, and Eq. I|3.1|l gives in both cases a parametrically larger rate 
w so ft(rn S oft/A) 2 for multiple scattering off of the non-perturbative background. 
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FIG. 4: A particle with momentum k S> m so ft scattering off of the non-perturbative field modes 
with wavenumber m so f t . 
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FIG. 5: 2— >1 processes catalyzed by a background field in (a) scalar ^ 4 theory and (b) gauge theory 
(inverse bremsstrahlung) . 



background of an oscillating inflation field [12l|2l|. However, processes analogous to Fig. 0] 
were not relevant in that case because the background field carried zero momentum and so 
could not directly transfer momentum to the particles. In that case, the dominant process 
was Fig. Et, which is not as efficient, giving a smaller value of v in those applications. For 
our application, we consider in the Appendix the related process of Fig. |SJd: 2— >1 inverse 
bremsstrahlung of particles with momenta k 3> m so f t , catalyzed by the non-perturbative 
background field. We find that, due to the Landau- Pomeranchuk-Migdal (LPM) effect, this 
process is not competitive with multiple scattering via Fig. HJ 

Failing to find a more efficient process than multiple scattering off of the non-perturbative 
background, our provisional conclusion is that the spectral index v should be 2. This is 
consistent with results from simulations 0], which extract v = 2 ±0.2 from Coulomb-gauge 
spectra of the cascade. From (|2.7|) . we also get that the UV end of the cascade should grow 
as 



A(t) oc *V(4-) „ 

Using (jl.2|) . this means that our classical treatment (fy ^> 
t < l/g 2 m so{t . 



(3.4) 

1 for k <C A) is valid for times 



IV. SCREENING, SCATTERING, AND INTERNAL CONSISTENCY 

Having proposed a cascade with 



(4.1) 
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we should now ask whether the development of the cascade might have a back-reaction on 
the dynamics of the infra-red (IR) modes with wave- numbers of order m so f t . These IR modes 
are the starting point for energy to enter the cascade, and one could imagine that a back- 
reaction might conceivably change the rate at which energy is taken into the cascade from 
the original anisotropic hard particles. If such a change to the rate depended on A(t) and 
so on time, it would be harder to understand why the energy in the cascade grows linearly 
with time, as indicated by numerical simulations. 

Let us investigate how the cascade affects the soft scale m so f t itself, given by (jl.3J) . This 
scale is initially created by the hard particle distribution, as we have indicated explicitly in 
(jl.3j) . However, once the cascade develops, its excitations could also contribute: 



where k is integrated over momenta in the cascade. How big is this contribution, and is 
its dependence on time through A(t) significant? For v < 2, the integral (J4.2|) would be 
sensitive to large A(t); for v > 2, it would not. The spectral index v = 2 found in this paper 
is the borderline case, and (j4.1|) produces 



If we ignore the logarithm, then Am 2 ~ ^soft> anc ^ the picture can be self-consistent. The 
creation of the cascade would make an 0(1) relative change to the IR physics driving the 
instability, and so could make an 0(1) relative change to the rate of linear energy growth. 
That still allows the system to settle down into steady growth as t — > oo, consistent with 
numerical results. 

The formal logarithmic dependence of ()4.3|) on A(t), however, suggests that the energy 
growth rate should see logarithmic in time corrections. Current simulations [18] do not 
support such logarithmic corrections, perhaps because the coefficient happens to be small, 
and perhaps because of the limited statistical power of existing simulations. In addition, 
the above analysis was rather crude: the details of exactly how Am 2 affects instability 
development depend, for instance, on the degree of anisotropy of the cascade particles as 
a function of time. 4 Thoroughly understanding the physics of the cascade at the level of 
logarithms remains an open problem. 



4 For instance, isotropic distributions do not produce magnetic instabilities. Adding a large isotropic con- 
tribution j' cascadc to an anisotropic hard particle distribution J hard will not change the momentum scale of 
the instability. It will, however, decrease the perturbative growth rate of such instabilities. (See, for exam- 
ple, the discussion in Rcf. 4].) Roughly speaking, this is because changing magnetic fields require electric 
fields, the ability of a medium to support electric fields depends on its conductivity, and the conductivity 
is affected by the number of isotropic as well as anisotropic particles. On the other hand, if J cascado is 
significantly anisotropic, it could increase the momentum scale of the instability, which would tend to 
increase growth rates. Numerical simulations sug gest that the excitations in the cascade are somewhat 
but not perfectly isotropic. A simulation in Ref. 18], for a particular hard particle distribution, found 
that the magnetic fields making up the excitations in the cascade had B^ equal to roughly 80% of B^. and 




(4.2) 




(4.3) 
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We see that cascade particles can have an 0(1) effect on the infrared dynamics (up 
to logarithms). In previous work j24[, we argued that particles can make two important 
contributions to the behavior of the plasma. One is the effect on infrared fields ( "screening," 
in the parlance of Ref. [24]). The other is the contribution to the scattering rate for hard 
particles. That is, when a hard particle (or another cascade particle) scatters, what fraction 
of the time is it scattering from a cascade particle? This is determined by the relative 
contributions to the integral 



Jk 

where (1 + /&) is a final state Bose stimulation factor and where we will focus on fc's in the 
cascade. For our / oc k~ 2 cascade (with / ^> 1), the integral is of form 



which is dominated by the smallest momenta and insensitive to the moving cutoff A(t). That 
is, the cascade particles are not important as scatterers. Most scattering is off the m soft 
nonperturbative background, not off the higher-momentum fluctuations. This is another 
way of seeing why it is scatterings from the soft background, rather than between cascade 
particles, which determine the progress of the cascade. 

V. WHITHER BOTTOM-UP THERMALIZATION? 

We introduced this paper with a discussion of how nonabelian plasma instabilities play a 
significant role in the bottom-up scenario for quark-gluon plasma equilibration in the high 
energy limit. We now have a reasonable picture of the development of non-Abelian plasma 
instabilities in a non-expanding system. Once we understand, at least parametrically, the 
basic processes and their rates, we should be able to apply this knowledge to an analysis 
of expanding plasmas. Unfortunately, there is an obstacle. For reasons of computational 
practicality, simulations have been performed for moderately, but not extremely, anisotropic 
hard particle distributions. For moderate anisotropy, there is a single parametric scale m so ft 
for IR physics. Non-perturbative IR magnetic fields are therefore 0(m 2 o{t /g 2 ), the linear 
growth of soft energy density is O(m^ oft t/g 2 ), and so forth. However, as we shall review, 
the first stages of the bottom-up scenario involve extremely anisotropic hard particle distri- 
butions, where (in local fluid frames) the angle A8 between typical particle velocities and 
the transverse plane is parametrically small. One must therefore understand the parametric 
dependence of instability development on this additional small parameter, A9. We leave 
such understanding for future work. In this section, we simply wish to emphasize the need 
by showing how the current state of uncertainty can have a significant impact on the first 
stage of bottom-up thermalization. In particular, we will consider a possible alternative to 
the recent analysis of the effects of instabilities by Bodeker 0]. 

In their original analysis of the stages of thermalization for arbitrarily high energy heavy- 
ion collisions, Baier et al. jl| considered an initial, non-perturbatively dense state of saturated 
low-x gluons produced in the collision. In their analysis, these were labeled the "hard" 
partons in the plasma, with momenta phard of order the saturation scale Q s . Their analysis 
is restricted to following the process of thermalization for r ^> 1/Q S , where r is proper time 
since the collision. 




(4.4) 




(4.5) 
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They assume initial conditions at time r ~ 1/Q S with number density n ~ Q^/ 1 g 2 of 
such hard particles. This number density dilutes as n ~ {Q s r)~ l Ql/ 1 g 2 due to the one 
dimensional expansion of the system. The scale (|1.3|) characterizing infrared physics (which 
they interpreted as a Debye screening scale) then evolves with time as 

m - ~ ? \ k ~ 7T ~ WW 1 ■ ( 51 ) 

Jp P 

Our notation u m 00 " for this scale is in preparation for our later discussion of plasma insta- 
bilities, for notational consistency with previous work. 5 

In the first phase of their analysis, the gluons dilute to perturbative densities as the system 
expands, and their local distributions of velocity temporarily become highly anisotropic due 
to momentum selection: over time, approximately free-streaming particles with the same v z 
would move to approximately the same displacement z, along the beam axis, from the origin 
of the collision. This drive towards local anisotropy is moderated only by processes which 
scatter the hard particles and so try to isotropize them. Baier et al. considered random, 
individual, small-angle 2—>2 scattering of the hard particles and found that, early in their 
scenario, the balance between expansion and scattering effects produces local distributions 
with anisotropy 

K T&rW (5 - 2) 

in local fluid rest frames. This result indicates extreme anisotropy for Q s t 3> 1. Later, 
however, once the system dilutes enough that /hard < 1, Baier et al. found that other 
processes come into play which eventually locally isotropize and thermalize the system. 
Here, we modestly focus our attention on the first phase, described by (J5.2j) in the original 
scenario. 6 

Arnold, Lenaghan, and Moore pointed out that the first phase will be modified by 
plasma instabilities. Recently, Bodeker 0] has attempted to deduce how scattering of hard 
particles from soft fields created by plasma instabilities will modify (|5.2|) . Bodeker assumed 
that soft fields consist only of a non-perturbative component (/soft ~ l/^ 2 ) with momenta 
?7tsoft ~ f^oo given by (jl.3j) . For hard particle distributions with 0(1) anisotropy, this is 
parametrically the only characteristic scale of unstable modes. With these assumptions, 
Bodeker found 

Pz_ (^SQft^) 1/2 1 /.ON 

P±~ Q s ~(Q.t)V*' 1 1 

which is parametrically larger than (|5.2|) . This estimate is in keeping with the discussion of 
this paper; the dominant way in which high momentum particles (whether in the cascade, 
as we have discussed, or even higher momentum, as is relevant here) undergo scattering, is 
by interaction with the nonperturbative fields at the scale m so ft- 

However, this estimate is not secure, because the treatment of the unstable scale m so f t 
has been made assuming 0(1) anisotropy, yet it predicts a parametrically large anisotropy 
for times r ^> 1/Q S . Therefore we must reconsider the behavior of the instability for highly 
anisotropic systems, a point also emphasized by Bodeker 0]. As we shall review, there 



5 The historical reason for the subscript "oo" is that, in addition to its role in the physics of instabilities, 
moo turns out to be the mass scale of transverse plasmons in the large momentum limit (i.e. p — > oo 

) 0- 



—3 /2 

° This corresponds to 1 <C Q s t a s ' in the original scenario of Ref. 
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is another characteristic scale g max ^> moo for such extremely anisotropic distributions. 
We believe that Bodeker's estimate of taking m so f t of order the lower scale (moo) in ()5.3|) 
probably provides a lower bound on the result for p z /p. We will not attempt to make a 
complete treatment here, but will provide what we believe is a sensible parametric upper 
bound. 

For highly anisotropic, oblate distributions of hard particles, the spread of angles of hard 
particle momenta p about the transverse plane is characterized by 

A# ~ — . (5.4) 
V 

The perturbatively unstable modes are described by two scales. One is moo, defined before, 
whose inverse represents the time scale that a gauge field must act on the particles before 
the back-reaction of particles on the field becomes important, moo characterizes the growth 
rate of instabilities. The other scale is the largest wave vector q of any unstable soft gauge 
field mode, which is 

gmax ~ ■ (5.5) 



For a detailed explanation of how these scales arise, see Ref. |2j; here we give a physical 
argument. For a gauge field mode to be unstable, the hard particles must fly for a time 
scale ~ in the presence of that mode, and be deflected in the same direction that whole 
time. That means that the magnetic field must be roughly the same all along the trajectory 
of the particle's path. Therefore, the gauge field must take roughly the same value over 
distances 5x± ~ 1/moo in the two transverse directions, and over a distance 5z ~ AO/rrioo in 
the longitudinal direction (since this is how far a typical particle moves in the longitudinal 
direction in that time scale). That means that modes with q± < moo and q z < vn^j 'A8 are 
unstable. This agrees with a more detailed analysis j2[. 

How large can the soft fields grow? We can place an upper bound as follows. Gauge fields 
in different color directions are all growing simultaneously. When these grow large enough 
that they can randomize the color of a hard excitation in a time scale ~ 1/moo, then the 
hard modes will no longer be contributing to the time development of the unstable modes 
on the time scale required to make the instability work. The rotation of color charge as a 
hard particle moves along a trajectory is a matrix given by the Wilson line 



"Pexp 



ig / dx ■ A 



(5.6) 



where V represents path ordering. In the analysis of a particular unstable mode, think 
of all the other unstable modes as background. The color rotation (|5.6|) due to a non- 
commuting background field will therefore be significant, and so stop instability growth, 
if g 5x ■ A > 1. To avoid color randomization of generic hard particles in our extremely 
anisotropic distribution, we therefore need, for example, 7 ' 8 

A± < ~ (5.7) 

gdx± g 



7 In order to discuss parametric estimates concerning the gauge field A, one should assume that one is in a 
gauge where A is relatively smooth. Here, we imagine that we study the state of the system at a particular 
time t by (i) working in Aq = gauge, and then (ii) making a spatial gauge transformation to make A as 
smooth as possible at that particular time. 

8 Another way to arrive at the relation l|5.7|l is to consider the propagator (D + v ■ D)^ 1 of particle 
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Since the dominant instabilities have q nearly in the z direction, with |q| ~ oWx, this implies 
that an upper bound on how large the corresponding magnetic field strength can grow is 9 

2 

coherent on transverse length scales of m^ 1 . Such fields will give momentum kicks to the 
hard modes due to the gv x B term in the Lorentz force law. The force is coherent over the 
time scale 5t ~ m^ 1 (since this is both how often B changes as one moves in the transverse 
direction, and how often the particle's charge gets rotated). The individual momentum kicks 
are of size Ap z ~ gB5t ~ q max , and they lead to momentum diffusion in the p z direction of 
size 

vl ~ j t (Ap*) 2 ~ gLx^ooT. (5.9) 

Therefore, 

l>. ^ Q, ~ QsA6 ■ ■ ' 

[Compare to (|5.3|) .] Substituting into the defining relation (|5.4|) for A9, and using the 
relation (|5.1|) for the time development of m^, gives 



EL^J^) 1 ' 2 L_ (5.11) 

The above estimate is correct if our picture of how plasma instabilities are cut off in a 
highly anisotropic setting proves correct. While we believe that our estimates represent an 
upper bound on such a cutoff, we are not confident that they are correct. The moral of this 
story is that better theoretical understanding is needed of the non-perturbative dynamics of 
instabilities in the case of extremely anisotropic distributions. Until then, it will be difficult 
to analyze even the first phase of bottom-up thermalization with complete certainty. 



VI. CONCLUSIONS 

We have argued that, in the presence of non-Abelian instabilities that have grown non- 
perturbatively large, the spectrum of the cascade of soft field energy from unstable modes 
into the ultra-violet should have the form ~ k~ u with v = 2. This is consistent with 
results from numerical simulations. 

We should emphasize that this turbulent cascade is only a transitory phenomenon ap- 
pearing during the local thermalization of the quark-gluon plasma in the high energy limit. 
It is not directly related to the spectrum of particles that much later leave the collision after 



fluctuations in non-Abelian Vlasov equations. In perturbative formulas for instabilities (e.g. Refs. [A lift 
11]), such propagators appear in the perturbative form (do + d- V)~ 1 ~ (ui — v ■ q)^ 1 . One can then ask 
when the gv ■ A in v ■ D can be treated as a perturbation to u> and v ■ q ~ q± + A9 q z , both of which are 
order m^. 

q± ~ vrioo. Similar reasoning to (|5.7|l could allow A z ~ l/(g5z) ~ q max /g, which gives an 0(q±A z ) 
contribution to B that is the same size as (15.811 . 
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hadronization. Also, though the word "turbulence" is used, the existence of this cascade 
does not mean that the plasma is described by ideal hydrodynamics at the time the cascade 
is formed. At the early times considered in this paper, most of the energy that will become 
the quark-gluon plasma is still in hard particles, Phard ~ Qs, which have not yet thermalized 
and which are locally extremely anisotropic. 

It would be delightful to be able to apply the current understanding of plasma instability 
development to produce a complete picture of thermalization in the high energy limit, fi- 
nally completing the program initiated by the original bottom-up scenario. A lot of progress 
has been made over the last year in understanding the qualitative behavior and paramet- 
ric dependence of the non-perturbative physics of instabilities for the case of moderately 
anisotropic hard particles. Unfortunately, figuring out the detailed role of instabilities in 
bottom-up thermalization will likely require at least a similar level of understanding of the 
case of extremely anisotropic hard particles, which have not yet been investigated non- 
perturbatively. 

Acknowledgments 

We would like to thank Larry Yaffe and Berndt Mueller for useful conversations. This 
work was supported, in part, by the U.S. Department of Energy under Grant Nos. DE- 
FG02-97ER41027, by the National Sciences and Engineering Research Council of Canada, 
and by le Fonds Nature et Technologies du Quebec. 

APPENDIX A: INVERSE BREMSSTRAHLUNG 

In this appendix, we consider the 2— >1 process of inverse bremsstrahlung, catalyzed by 
the soft, non-perturbative field, as depicted in Fig. We will treat a case slightly more 
general than that depicted in the figure, namely the case of one particle with momentum 
k ~ A and another particle with a momentum k! larger than m so f t , up to and including the 
scale A. We will find that the coherence time of this process is longer than the mean free 
time between collisions of the particles with the soft background (via Fig. HJ), and so we will 
have to account for the Landau-Pomeranchuk Migdal (LPM) effect [25( in estimating the 
rate of inverse bremsstrahlung. 

For an initial estimate, however, let us momentarily ignore the LPM effect and con- 
sider the diagram of Fig. El in isolation. Bremsstrahlung in a soft background is a small- 
angle, nearly-collinear process. Parametrically, the rate for bremsstrahlung or inverse 
bremsstrahlung may be estimated as (i) the rate for small-angle scattering from the back- 
ground field (Fig. 0J), times (ii) a factor of g 2 for absorbing or emitting the additional gluon 
in Fig. times (iii) an initial or final state factor of / or 1 + / for that gluon, and (iv) a mo- 
mentum integral dk! j k! (responsible for the logarithmically large rate of soft bremsstrahlung 
emission in vacuum). As discussed in Sec. IIII[ the rate for small-angle scattering from the 
non-perturbative background is O(m so ft)- This gives 



for (inverse) bremsstrahlung. To find the rate for a particle with momentum k ~ A to 
substantially change its energy, we must correct for two things. First, absorbing a momentum 




(Al) 
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k' < k makes only an order k'/k relative change to k. Second, it is the difference between 
the rates of absorption and emission of a particle of momentum k! which is relevant. This 
involves the difference f(k) — f(k + k') ~ f(k) k'/k. Therefore we must introduce two powers 
of (k'/k) into the estimate above, to get 

2 fdk'fk'V f A fm so{t y k'dk' 

for doubling energy through inverse bremsstrahlung. For v < 2, this is dominated by k' ~ A, 
and gives T ~ ml^A~ u (up to logarithms when v = 2). If we equated the power of A in 
(|A2|) with that in (|2.7|) . we would obtain z/ = 2, which would imply that this process was 
parametrically just as important as the multiple soft scattering treated in Sec. IIHI However, 
the rate (|A1J) is an overestimate; we must account for the LPM effect. 

In Abelian theories, LPM suppression occurs because (i) when a charged particle scat- 
ters, the probability that a photon will be bremsstrahlung emitted is not very sensitive to 
the scattering angle, and (ii) an emitted photon of a given wavelength cannot resolve the 
difference between a single collision and A" collisions if the A" collisions are close enough to- 
gether. 10 The largest such A" can be determined parametrically by equating (i) the time Nt 
for A" consecutive scatterings of the form of Fig. HJ where r is the mean free time between 
scatterings, with (ii) the formation time 1 / (5E) for the complete 1 —>■ 2 bremsstrahlung (or 
2 — >• 1 inverse bremsstrahlung) process, where 8E is the off-shellness of the energy during 
the process, which is 

6E=\E 1 + E 2 -E*\ (A3) 

for a 2<->l process with kik 2 <->• fc*. For the reasons of resolution already mentioned, 
each A" collisions are equivalent to a single collision in terms of the probability to induce 
bremsstrahlung, and so the rate Y for (inverse) bremsstrahlung processes is a factor of A" 
smaller than the naive estimate (jAl|) which treats every scattering as independent. 

The same sort of estimate works for gluon bremsstrahlung in non- Abelian gauge theories. 
Since bremsstrahlung is an almost collinear process, we can take k± <C k, where _L is relative 
to the collinear axis, and approximate (|A3|) for k ^> m so f t as 



2ki 2k 2 2K 

The main difference between photon bremsstrahlung in QED and gluon bremsstrahlung 
in QCD is that, in QCD, any of the three particles (including the bremsstrahlung gluon) 
can undergo the scattering. Therefore we may take each k\ ~ Nm^ oft , in which case SE is 
dominated by the smallest energy particle, which is k' in our case. Further, m 2 ~ <C k\ 
can be dropped, giving 

5E ~ S ~ . (A5) 

k k 

Plugging into ()A4|i . and equating the scattering time 1/SE to Nt as outlined previously, we 
have 

1 k' 

Nt ~ — ~ Ar 9 . (A6) 
5E Nml h 



10 For a review of the LPM effect in photon bremsstrahlung, see Ref. For a discussion by the present 

authors of LPM suppression in non- Abelian plasmas near equilibrium, see Refs. [j^. I27 1 . For a partial 
selection of earlier discussion of the LPM effect in non- Abelian gauge theories, see Refs. 
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Solving for N, and using the fact that r ~ l/m so f t for the individual scattering processes of 
Fig. HI determines 

N ~ Vfc'/msoft. (A7) 

Therefore, the rate is reduced by a factor of {m^/k') 1 ^ 2 ^ which is a large suppression for 
k' ^> m so ft. This modifies Eq. (jA2|) to, 

r fc^2fc,(LPM) ~ m soft / [-jr ) -p- • (A8) 

This expression is dominated by small k' if v > 3/2. The k' ~ A edge of the range of 
integration provides T ~ A - " -1 / 2 , which combined with (|2.7|) would now yield f = 7/4. 
Therefore, hard bremsstrahlung is not a competitive process. However, the k' ~ m so ft edge 
of the integration range gives T ~ mg oft /A 2 , which combines with Eq. (|2.7|) to give v = 2. 
Absorption of an excitation with energy of order m so f t is competitive with (and possibly 
indistinguishable from) scattering from such excitations, but inverse bremsstrahlung of all 
harder scales is sub dominant. 
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